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THE  THERMAL  REGIME  OF  THE  DOWNSTREAM  SHOULDER  OF  ROCKFILL  DAMS 


A ays tea  of  equations  describing  the  thermal  regime  of 
rockf ill  dams,  one  taking  natural  air  convection  into 
account,  is  presented  in  the  paper.  A numerical  solu- 
tion of  the  equations  is  given  for  the  profile  of  the 
downstream  shoulder  of  the  rockflll  dam  of  the  Vilyuy 
Hydroelectric  Power  Plant,  the  solution  being  obtained 
with  the  aid  of  a computer.  The  results  of  calculation 
and  of  field  observations  are  found  to  be  in  close  agree- 


Thls  paper  is  devoted  to  study  of  the  thermal  regime  of  the  downstream 
shoulder  of  rockflll  dams.  Calculation  of  the  thermal  regime  of  the  down- 
stream shoulder  is  necessary  in  order  to  establish  the  thermal  state  of  the 
apron  (core),  the  base  of  the  dam,  and  the  location  of  drainage  Installations. 
The  thermal  regime  of  the  downstream  shoulder  of  a rockflll  dam  is  formed  as 
a result  of  conductive  heat  transfer  and  convective  heat  transfer  by  the  air 
moving  in  the  pores. 

In  order  for  movement  of  air  to  occur  in  the  body  of  the  dam  it  is  neces- 
sary for  the  lifting  force  (due  to  the  heterogeneity  of  the  temperature  field) 
to  be  greater  than  or  equal  to  the  force  of  resistance  of  the  rock  fill  to 
the  movement  of  the  air.  The  rock  fill  causes  losses  due  to  deformation  of 
the  moving  air,  i.e.,  losses  from  contraction  and  expansion  of  the  air  stream 
and  losses  due  to  air  friction. 


The  ntaiber  of  contractions  and  expansions  per  unit  length  of  the  rock 
fill  depends  on  the  dimensions  of  the  jointings  and  their  geometric  shape, 
i.e..  Is  a function  of  the  diameters  of  the  jointings  and  of  the  porosity 
of  the  fill. 


The  force  of  the  resistance  due  to  friction  is  the  most  complex  in  na- 
ture. This  Is  to  be  ascribed  not  only  to  the  roughness  of  the  surface  of 
the  rocks  and  their  varying  orientation  relative  to  the  flow,  but  also  to 
the  temperature  of  the  moving  air  Itself  and  of  the  surface.  Hence  in  order 
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to  determine  the  force  of  friction  It  1*  necessary  In  the  general  case  to 
solve  together  the  problem  of  the  thermal  state  and  that  of  movement  of  the 
air,  making  allowance  for  the  boundary  conditions  on  the  surface  of  the  joint. 

In  practice  It  Is  not  possible  to  solve  this  problem,  especially  for  a 
rock  fill  and  its  variegated  structure.  In  view  of  this  circumstance  it  Is 
advisable  to  introduce  quantities  characterizing  the  phenomenon  on  the  average. 
For  example,  the  force  of  resistance  of  a rock  fill  is  regarded  as  a generalized 
force  which  can  be  determined  empirically. 

Transfer  of  heat  to  the  jointings  of  the  fill  is  effected  in  various 
ways.  At  the  places  of  direct  contact  between  two  adjacent  jointings  and 
through  interlayers  filled  with  stationary  air,  heat  is  transferred  by  thermal 
conductivity,  and  by  heat  exchange  in  the  areas  swept  by  the  moving  air. 

Hence  the  teaperature  even  of  a single  jointing  varies  from  point  to  point 
of  the  latter.  In  view  of  this  fact,  the  variation  in  the  heat  content  of  a 
rock  fill  at  any  point  may  be  expressed  through  the  variation  in  a certain 
mean  temperature  determined  from  the  thermal  balance  conditions. 

The  temperature  in  the  core  of  the  moving  flow  (with  the  exception  of 
the  wall  boundary  layer)  may  in  practice  be  assumed  to  be  the  same  over  the 
entire  cross-section,  since  the  air  moves  in  porous  channels  of  infinitely 
varying  cross-section  creating  the  conditions  for  mechanical  displacement. 

Hence  a flow  of  air  depending  on  the  thermal  regime  of  the  dam  itself 
is  formed  in  the  downstream  shoulder  of  a rockf ill  dam.  Consequently,  a flow 
of  air  variable  in  time  will  be  observed  until  the  temperature  of  the  down- 
stream shoulder  of  the  rockf ill  dam  is  everywhere  near  the  temperature  of 
the  external  air. 

The  known  methods  of  calculating  the  thermal  regime  of  the  downstream 
shoulder  of  rockfill  are  based  on  utilization  of  reduced  (effective) 

coefficients  of  thermal  and  temperature  conductivity,  which  were  derived 
earlier  for  one-dimensional  problems  [1],  [2],  [3],  i.e.,  for  a condition 
such  that  the  vectors  of  the  thermal  flow  and  of  the  filtration  speed  of 
the  air  in  the  pores  coincide  in  the  region  under  consideration. 

The  geometric  dimensions  of  the  downstream  shoulder  of  any  rockfill  dams 
are  found  to  be  such  that  a moving  flow  of  air  cannot  be  described  as  one- 
dimensional,  but  rather  must  be  described  in  at  least  two  spatial  coordinates 
(x,y).  Owing  to  the  two-dimensionality  of  the  flow  in  the  pores  of  the  fill, 
the  thermal  regime  of  the  higher  lying  jointings  will  be  determined  by  thermal 
calculation  of  all  the  lower  lying  layers,  and  vice  versa. 

As  is  demonstrated  by  experience  acquired  in  erecting  rockfill  dams,  the 
grain-size  distribution  of  a rock  fill  varies  over  a wide  range.  The  larger 
pores  of  a fill  usually  are  filled  with  smaller  fractions,  i.e.,  the  area  of 
the  particles  filling  the  unit  of  volume  of  the  fill  is  a large  one.  Conse- 
quently, the  forces  of  friction  will  be  larger  than  the  forces  of  inertia  of 
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the  moving  air.  Hence  the  flow  of  air  In  the  pores  may  be  regarded  as  a 
filtration  process  [4]. 

1.  Mathematical  Statement  of  the  Problem 

If  It  is  assumed  that  the  rock  fill  is  isotropic,  the  porosity  does  not 
vary  with  time,  phase  shifts  are  absent  from  the  water  vapor  forming  part  of 
the  components  of  the  moving  air,  the  conductive  thermal  conductivity  of  the 
air  in  the  direction  of  movement  is  small  in  comparison  to  the  convective 
heat  transfer  and  may  be  disregarded,  the  volumetric  heat  capacity  of  the 
rock  fill  is  constant  and  does  not  depend  on  the  coordinates,  other  forma 
of  mass  transfer  (such  as  water  filtration  or  infiltration)  are  absent,  and 
the  variation  in  air  density  in  time  and  the  vertical  variation  in  barometric 
pressure  are  small  [5],  then  the  following  relations  are  valid  for  the  moving 
air  per  unit  volume: 

(1) 

(2) 

(3) 


(4) 

where  -p-  are  the  pressure  gradients  respectively  along  the  ox  and  the 
y oy  axis ; 

F(u),  F(v)  are  the  pressure  drops  due  to  air  friction  in  the  fill  re- 
spectively along  the  ox  and  the  oy  axis; 

m is  the  porosity  of  the  fill; 

6 is  the  coefficient  of  volumetric  expansion  of  air; 
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y is  the  volumetric  weight  of  air; 

AT  ■ 0 - t is  the  temperature  difference  of  the  pore  and  the  external 
ext  air; 

p is  air  density; 

Pq  is  air  density  at  zero  temperature; 

6 is  the  temperature  of  the  pore  air; 
x,  y are  coordinates; 

u,  v are  the  speed  of  air  respectively  along  the  ox  and  the  oy  axis; 
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t is  ths  temperature  of  Che  external  air. 

Let  us  reduce  equations  (1)  and  (2)  to  a form  convenient  for  calcula- 
tion. For  this  purpose  we  differentiate  equation  (1)  with  respect  to  y and 
equation  (2)  with  respect  to  x. 

Then  we  obtain 

d'P  dF(u) 
dxdy  dy 

d'P  dFlv)  d(m,37Ar) 

dx  dy  dx  dx 

or  dF(u)  _ dfc(v)  _ d(mfa±T) 

dy  dx  dx 

A larger  number  of  papers  have  been  devoted  to  the  problem  of  determina- 
tion of  the  pressure  drops  in  porous  media  on  movement  of  gases  and  liquids 
through  them.  Many  authors  regard  flow  through  a porous  layer  as  an  inter- 
nal problem  (flow  in  channels  between  particles  [6],  [7]),  and  others  as  an 
external  one  (flow  of  air  around  a sphere  [8]). 

On  the  basis  of  these  papers  various  empirical  and  semiempirical  formu- 
las have  been  proposed,  ones  generalizing  a large  number  of  experiments  con- 
ducted with  particles  of  small  diameter  in  one-dimensional  flow.  The  rate 
of  filtration  enters  these  relations  in  the  second  power;  this  is  not  ac- 
ceptable for  describing  losses  in  a two-dimensional  region,  since  a linear 
relationship  must  exist  between  the  velocity  vector  and  its  components  in 
a two-dimensional  region  [9]. 

Hence  it  is  necessary  to  adopt  the  general  expression  of  the  Darcy 
law  [10]  to  describe  the  losses  due  to  friction: 

F(u)  “ -k  (pu)  and  F(v)  “ -k  (pv),  (8) 


(5) 

(6) 
(7) 


where  k - v/k'; 

v is  the  coefficient  of  kinematic  viscosity  of  air; 

k'  is  the  coefficient  of  permeability  of  the  fill. 

Inserting  (8)  into  (7) , we  obtain 

d(kpu)  d(kpv)  _ d(mrft±T)  (9) 

Qy  ox  ox 

Equation  (9)  must  be  considered  in  conjunction  with  equation  of  dis- 
continuity (3).  It  is  convenient  for  practical  calculations  to  reduce 
equation  system  (9)  and  (3)  to  a single  equation. 

For  this  purpose  it  is  necessary  to  express  the  mass  velocities 
through  the  flow  line  function: 


Equation  (9)  then  assumes  the  form 


0 (/nplA  T) 
Ox 


This  Is  an  elliptical  equation  with  variable  coefficients.  Equation 
(11)  la  also  valid  for  an  anisotropic  medium  If  the  axes  of  the  anisotropy 
coincide  with  the  coordinate  axes,  but  with  the  porosity  constant.  The 
equation  of  discontinuity  may  subsequently  be  disregarded,  since  It  is 
converted  to  identity. 


Thus  we  obtain  the  following  system  of  equations  with  boundary  con' 
ditions  (Figure  1)  with  which  to  determine  the  thermal  regime  of  the 
downstream  shoulder  of  a rockf ill  dam: 


(*?-rwrC'  a 
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Figure  1.  Illustrating  determination  of  boundary  conditions:  1.  UVB 

2.  Core 


(1)  the  temperature  on  the  day  of  the  downstream  shoulder 


q when  * « *ur  F = W 


(2)  on  the  boundaries  of  conjugation  of  the  downstream  shoulder  with 
other  airtight  materials 
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(3)  the  temperature  of  the  moving  air  on  the  boundaries  of  the  down- 
stream shoulder  are: 

(a)  eoc,  y,  T)-e<*sur,ysuI^whM  x - xsor>  y - ysur> 

(b)  6(x,  y,  T)  -9(xr.s.yre/)when  . - xrei>,  y - 

(A)  the  flow  line  function  is 

(a)  m 0,  when  0 ^ x £ B,  y ■ 0 and  0 < y £ H,  x - 0, 

(b)  ^ (xi  y,  t)  « i|>  (x  y ,t)when  x - x , y - y 

v/  r ' 7 r xur,  sur  sur  7 7sur 

With  the  Initial  conditions 

t (x,  y,  *)  = <(*.  y,  0); 

6 (a-,  y,  v)  =>  0(-t,  y,  0); 

(jc,  y,  t)  = •!>  (*,  y,  0). 

Here  t Is  the  mean  temperature  of  the  rock  fill; 

tsur  is  the  temperature  on  the  surface  of  the  rock  fill; 

9aur  is  the  temperature  of  the  moving  air  on  the  surface  of  the  dam; 

T is  time; 

is  the  coefficient  of  volumetric  heat  exchange; 

CvQ^  is  the  volumetric  heat  capacity  of  the  rock  fill; 

C . is  the  volumetric  heat  capacity  of  the  air; 

VOX 

q is  the  amount  of  heat  transferred  to  the  interior  of  the  dam  by 
thermal  conductivity; 

If 

X -r—  Is  the  thermal  flux  from  the  direction  of  the  downstream  shoulder 
n toward  the  line  of  conjugation; 

3t 

X ,■ ' re8  is  the  thermal  flux  toward  the  line  of  conjugation  from  the 

res  n direction  of  the  base  and  core  (apron); 

X is  the  coefficient  of  thermal  conductivity  of  the  rock  fill  without 
convection; 

x , y are  the  coordinates  of  conjugation  of  the  downstream  shoulder 
res  res  tbe  base  and  core; 

Xrefl  is  the  coefficient  of  thermal  conductivity  of  the  base  and  apron; 


(18) 

(19) 

(20) 

(21) 
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In  order  to  determine  9(xgur>  ysur»T^  and  ^^xSur'  y8ur'  on  the  day» 

together  with  the  equations  of  system  I It  Is  necessary  to  consider  equations 
describing  the  velocity  and  temperature  fields  of  the  air  outside  the  dam. 


It  is  a characteristic  of  the  process  of  natural  convection  in  the 
downstream  shoulder  of  a rockf ill  dam  that  the  warm  air  emerging  is  pressed 
against  the  core  (apron)  of  the  dam.  This  is  explained  by  the  fact  that 
the  core  (apron)  of  the  dam  has  a higher  temperature  than  do  the  other  parts 
of  the  downstream  shoulder.  Hence  escape  of  air  always  takes  place  on  the 
crest  of  the  dam  (if  there  is  no  obstacle  inside  the  dam  in  the  form  of 
airtight  interlayers) . The  warm  air  escaping  from  the  dam  has  no  appreci- 
able effect  on  the  temperature  regime  of  the  surrounding  external  air  of 
the  downstream  shoulder  of  the  dam,  as  has  been  established  by  observations 
at  the  Vilyuy  Hydroelectric  Power  Plant. 
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Figure  2.  Arrangement  of  thermistors  on 
the  downstream  slope  of  the  dam  of  the 
Vilyuy  Hydroelectric  Power  Plant 


The  results  of  the  observations  are  illustrated  in  Figure  3.  As  is 
to  be  seen  from  Figure  3,  similar  variations  take  place  in  the  air  tempera- 
ture on  the  surface  of  the  dam.  However,  at  point  No.  2 (Figure  2)  there 
is  observed  an  appreciable  movement  of  air  into  the  interior  of  the  dam, 
the  temperature  of  which  equals  the  external  air  temperature. 


At  a distance  of  1.0-1. 5 m from  the  point  of  escape,  the  temperature 
of  the  air  moving  in  the  pores  is  near  the  external  air  temperature  (Figure 
3) ; the  variation  in  the  former  temperature  does  not  exceed  4s . Hence  the 
value  of  3 (mByAT)  on  the  day  may  be  determined  from  the  difference  between 
3x 

the  air  temperature  on  the  surface  of  the  downstream  shoulder  and  the  ex- 
ternal air  temperature.  An  error  has  an  effect  at  a short  distance  from 
the  point  of  escape  of  the  warm  air  from  the  dam,  within  a range  of  5 m. 


Figure  3.  Temperature  of  external  air  on 
dam  surface  according  to  readings  of  mercury 
thermometers.  1.  T,  days.  2.  November, 
1968. 


H3s6fH  IS  St  l 


In  this  case  the  velocity  field  of  the  air  on  the  day  of  the  down- 
stream shoulder  of  the  dam  can  be  determined  from  the  physical  nature  of 
the  problem,  which  assumes  that  the  velocity  of  approach  (departure)  of 
the  external  air  equals  the  filtration  velocity  of  the  air  in  the  pores 
on  the  day  of  the  downstream  shoulder.  This  condition  is  fulfilled  if 
it  is  assumed  that  d\p/dn  ■ 0,  i.e.,  the  movement  of  the  air  is  always 
normal  to  the  surface  on  the  day  of  the  downstream  shoulder.  The  assump- 
tions adopted  permit  solution  of  the  problem  of  the  thermal  regime  of  the 
downstream  shoulder  of  a roclcflll  dam  without  analysis  of  the  velocity 
and  temperature  fields  of  the  air  outside  the  dam.  It  should  be  noted 
that  the  velocity  and  temperature  fields  of  the  air  outside  the  dam  de- 
pend on  many  meteorological  factors  which  we  have  not  as  yet  succeeded 
in  alloying  for. 

Equation  system  I is  also  valid  for  regions  having  a more  complex 
geometric  shape,  and  in  particular  for  the  downstream  shoulder  of  a rock- 
fill  dam  with  an  apron.  In  this  case  as  well  the  external  air  temperature 
i.e.,  AT  ■ 0 - t , must  be  taken  as  the  origin  of  the  temnerature  level 
[13]. 

The  results  obtained  were  employed  to  calculate  the  thermal  regime 
of  the  downstream  shoulder  of  the  rockflll  dam  of  the  Vilyuy  Hydroelectric 
Power  Plant. 

To  establish  the  initial  and  boundary  conditions  it  is  necessary  to 
know  the  climatic  and  cryopedologic  conditions  of  the  transit  line  of  con- 
struction of  the  dam,  the  filtration  properties  of  the  base,  the  coarse- 
ness of  the  jointings  and  the  porosity  of  the  rockflll,  and  the  structure 
of  the  dam  itself. 

2.  Climatic  and  Crvooedoloaic  Conditions  of  the  Area  of  Construction 


The  area  of  construction  of  the  Vilyuy  Hydroelectric  Power  Plant  is 
characterized  by  a long  cold  winter  with  little  snow  and  a hot  summer. 

The  amplitude  of  the  annual  temperature  variations  reaches  100°  C,  the 
minimum  temperature  in  winter  -60°C,  and  the  maximum  temperature  in  sum- 
mer +36“  C.  The  winter  is  up  to  seven  months  long.  The  mean  annual  long- 
term temperature  of  the  area  is  -8.2°  C.  The  thickness  of  the  snow  cover 
does  not  exceed  250-300  mm.  The  annual  distribution  of  the  external  air 
temperature  is  illustrated  in  Figure  4. 

The  area  of  construction  of  the  water  engineering  system  is  situated 
in  an  area  of  continuous  permafrost,  with  prevailing  temperatures  of  -3  — 
-5°  C.  The  thickness  of  the  permafrost  rocks  is  estimated  at  250-300  m. 

A talik  with  a positive  temperature  can  be  traced  in  the  bed  of  the  Vilyuy. 
The  width  of  the  strip  of  thawed  rocks  corresponds  to  the  width  of  the  river 
during  the  summer  low-water  period  [12]. 


Figure  4.  Graph  of  external  air 
temperature. 


The  original  slopes  of  the  Vllyuy  valley  differ  widely  from  each 
other  in  permafrost  conditions.  The  left  bank  is  made  up  of  perenially 
frozen  rocks  the  temperature  of  which  ranges  from  -4  to  -5*  C. 


The  rocks  of  the  right-bank  slope  are  characterized  by  low  temperatures 
reaching  -7*  C.  The  depth  of  seasonal  thawing  in  the  rock  ranges  from  3 to 
4 m [19]. 


3.  Geological  Engineering  Conditions  of  the  Construction  Site 


According  to  the  data  of  the  S.  Ya.  Zhuk  Gidroproyekt , the  bed  portion 
and  the  lower  slopes  of  the  Vllyuy  River  in  the  section  line  are  character- 
ized by  constant  geological  engineering  conditions.  The  normal  diabases 
are  widely  distributed  here.  The  best  preserved  and  most  monolithic  rocks 
occur  under  the  bed  and  the  portion  of  the  right  bank  adjoining  it. 


Jointing  is  associated  with  the  primitive  rock,  being  due  to  chilling 
of  the  intrusion  under  conditions  of  tectonic  stress.  The  fissures  divide 
the  diabaslc  massif  into  columnar,  sheet,  and  slab  jointings.  The  opening 
of  the  fissures  is  small,  from  fractions  of  a millimeter  to  1.0-1. 5 cm. 
From  the  viewpoint  of  jointing  the  rocks  are  classified  as  low-jointed, 
i.e.,  with  less  than  10  cracks  per  running  water.  The  coefficient  of 
fissure  cavitation  is  small,  failing  within  the  limits  of  0.5-1.0Z. 


According  to  the  data  of  experimental  pumpings,  the  coefficient  of 
filtration  of  the  base  is  estimated  at  0.1  to  2 m/day.  The  coefficient 
of  filtration  of  the  adjoining  left-bank  area  may  reach  10  m/day. 


The  base  of  the  Hum  is  represented  by  firm  diabases  with  a small  CO' 
efficient  of  filtration;  on  the  rlverbank  sections  they  are  covered  with 
friable  deposits  1.5-2  m thick.  The  friable  deposits  were  removed  when 
the  dam  was  erected  [19]. 


4.  The  Rockfull  Dam  of  the  Vllyuy  Hydroelectric  Power  Plant 


The  blind  rockfill  dam  has  a height  of  74  m and  a length  along  the 
crest  of  around  600  m [19].  The  horizontal  equivalent  of  the  slopes  varies 
in  height.  Berms  used  for  temporary  passages  and  permanent  roads  are  in- 
stalled at  Intervals  of  10-15  m on  the  slopes.  The  rubble  mound  of  the 


....  ■ • m 


dam  was  built  up  by  year-round  operation  by  the  pioneer  method  with  un- 
aorted  stone  obtained  in  the  exploitation  of  useful  cuts  or  special  quar- 
P18  transitional  zone  is  in  the  form  of  two-layers  of  reverseqfil- 

by^L«s0of^JecS^rockrc^8^ingf8  °f  8ra±n  8iZ6S  °'4°  and  °"l5°  ””  obtain®d 

equivalent  of  1:1.4,  .ud  the  up.tr™  slop,  one  of To  »rI-~ 
Idee  “Jd®5®inin8  tbe  8011  during  fluctuation  of  the  head  bay,  the  upstream 
®d*J  °J . the  ap5°n . I8  covered  ^th  a filter  above  the  apron,  of  a natural 

d?rec“on8ofVthe”h^n*  T°  ensure  8tatic  stability  of  the  apron  from  the 
irection  of  the  head  bay,  the  apron  is  weighted  down  with  a rock  fill. 

of  a liStS^?  ff  th!  apr°n  With  the  rock  base  18  effected  in  the  form 
tatioj  S the  COTeueJPlaCe  wlth  a cementation  culvert.  Cemen- 

reservoir  is  filled,  h*1™  ^ ***  fr0m  th®  culvert  is  schedule  before  the 

The  grain-size  composition  of  the  fill  is  sufficiently  homogenous 

^r«^r»^iS^2’}  sLh  sssHssz  «£«- 


tl«ht  are  dan8ar’  “d  ln  cases  even  air- 

Th?s  is  i 7ar  waa  foraed  on  PK  2~40  at  an  elevation  of  54.0  m. 

50  6 m £ J eSqap*  °f  warn  air  b®y°nd  the  berm  at  an  elevation  of 

50.6  m.  The  fill  is  anisotropic  within  the  limits  of  one  layer  i e the 

mean  diameters  and  porosity  of  the  fill  are  not  equal  in  the  upier'wd 

lower  parts  of  the  layer.  However,  the  mean  porosity  value  of  the  den  of 
the  Vilyuy  Hydroelectric  Power  Plant  la  in  our  opinion  approximate]*  30lf 

The  material  of  the  fill  possesses  the  following  Uhy si comechanic el 

3r025rs/^  !!iTtrlC  W*i8ht  approxi“tely  2.97  g/cm5,  specific  gravity, 
3.025  g/cm  , and  temporary  compressive  strength  on  the  average  1700 kg/ cm2. 

The  thermophysical  characteristics  of  the  fill  material  as  determined 
by  the  two-point  method  with  a cylindrical  sample  (R  • 0.032m  * • 0 055 
are  ®8  follow*:  coefficient  of  temperature  conductivity  (i)  0.004685 

m /h,  and  coefficient  of  thermal  conductivity  (l#k>  2.8  kcal/m-h-deg. 

. Tb®  values  given  for  the  parameters  will  be  used  in  what  follows  in 
determination  of  the  coefficients  appearing  in  the  calculations. 
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5.  System  of  Calculation  of  the  Thermal  Regina  of  the  Ro clef 111  Dam 
of  the  Vllyuy  Hydroelectric  Power  Plant 

The  following  system  of  calculation  was  adopted  to  determine  the  thermal 
regime  of  the  downstream  shoulder  of  the  rockflll  dam  of  the  Vllyuy  Hydro- 
electric Power  Plant. 

The  rock  fill  is  isotropic,  the  porosity  is  constant,  the  rock  base 
does  not  filter,  and  the  latter  is  frozen.  An  exception  Is  made  by  the 

upstream  side  of  the  apron  and  the  upstream  weighting,  where  thawed  soil 

with  a moisture  content  of  20Z  to  a depth  of  4.5  m was  assumed.  The  warm- 
ing effect  of  snow  on  the  slopes  of  the  downstream  shoulder  is  disregarded, 
the  snow  not  representing  an  obstacle  to  movement  of  air  into  the  Interior 
of  the  dam.  The  data  of  the  mean  monthly  air  temperatures  were  employed  in 
the  calculations . The  coefficient  of  volumetric  heat  exchange  is  assumed 
to  be  constant  and  Independent  of  the  coordinates  and  time.  The  profile 
of  the  dam  along  the  downstream  slope  is  assumed  to  be  straight,  in  con- 
trast to  the  actual  one,  which  is  in  the  form  of  a broken  line.  No 

analysis  is  made  of  the  velocity  and  temperature  fields  of  the  external 
air  outside  the  dam,  it  being  assumed  that  at  a short  distance  from  the 
surface  the  temperature  of  the  air  equals  the  temperature  of  stationary 
distribution  at  a great  distance  from  the  dam. 

It  was  assumed  that  the  variation  in  the  density  and  the  volumetric 
weight  of  the  air  in  the  pores  as  a function  of  the  coordinates  will  not 
play  a significant  role  in  the  velocity  distribution.  The  coefficient  of 
permeability  in  equation  (14)  is  taken  to  be  constant  and  to  characterize 
the  entire  downstream  shoulder  as  a whole,  since  the  granulometric  composi- 
tion of  the  fill  and  its  filtration  properties  cannot  be  definitely  indi- 
cated [12]. 

Equations  (13),  (14)  of  system  I then  assume  the  form 


! 


30  , 1 <?V  oB 1 dif  <?9  a,  „ _ g. 

‘ m ay  ax  m ax  ay  ' 

!f rn§g  00 

~0x*  "^"Sy3  k ~ ox ' 

where  g is  the  acceleration  of  gravity. 

Boundary  conditions: 

(1)  temperature  on  the  surface  of  the  downstream  shoulder: 


(22) 

(23) 


1 


Cvoi°aur  W + Cv0i  when 


X = x • 
suf 


sut 


(24) 


(2)  on  the  edges  of  conjugation  of  the  downstream  shoulder  with  the 
other  materials,  boundary  conditions  of  the  fourth  kind  must  be  fulfilled, 
1 . e . , 
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x **-  • 

res  dn  ’ 


t = t 


red 


(25) 


(3)  temperature  of  air  in  movement  on  the  boundaries  of  the  down- 
stream shoulder: 


а)  0(*.  y.  *)  — t(xre»  y res  T) 

when 

x - Xr.«-  y = yred 

б)  9(x,  y,  T)  = ® (-*8ur  ysur  ^ - 

* I 

when 

*-x.ut  y-y.ur 

if  the  velocity  vector  is  directed  into  the  interior  of  the  dam; 


(26) 


(4)  flow  line  functions: 


* = o, 

when 

X = X 

res 

y ~ y res 

'I*  =*0, 

when 

X 3 Jfji 

y-y*; 

-v  = o, 

when 

X =*  x,\ 

y = 

(27) 

II 

o 

when 

X — X 

sur 

y ~ ^sur 

where  xl,  yi  are  the  coordinates  of  the  horizontal  airtight  Interlayer 
Z at  an  elevation  of  54.0  m; 

x_,  y.  are  the  coordinates  of  the  horizontal  part  of  the  crest  of 
*■  the  downstream  shoulder  of  the  dam. 


The  initial  temperature  of  the  rock  fill  can  be  determined  on  the 
basis  of  the  data  of  field  observations  at  a specific  time,  i.e.,  t(x,  y, 
x)  - t (x,  y,  0). 

In  order  to  solve  system  I it  is  necessary  to  know  the  initial  tem- 
perature of  the  air  moving  in  the  pores.  No  measurement  apparatus  for 
determining  the  temperature  of  the  moving  air  was  installed  on  the  dam 
of  the  Vilyuy  Hydroelectric  Power  Plant.  Hence  the  initial  temperature 
of  the  moving  air  must  be  assumed,  with  a certain  amount  of  tolerance,  to 
equal  the  initial  temperature  of  the  rock  fill. 

Figure  5 illustrates  the  temperature  field  of  the  downstream  shoulder 
of  the  rockfill  dam  of  the  Vilyuy  Hydroelectric  Power  Plant  as  of  1 Novem- 
ber 1967,  this  being  taken  as  the  initial  condition. 

The  necessary  data  on  the  coefficient  of  heat  exchange  of  the  rock 
fill  with  the  air,  the  area  swept  by  the  air  stream  per  unit  fill,  the 
pressure  drop  due  to  resistance  of  the  rock  fill  to  the  movement  of  the 
air,  and  the  coefficients  of  thermal  and  temperature  conductivity  with 


I 
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do  allowance  for  convection  can  be  determined  on  the  basis  of  field  observa- 
tions and  laboratory  studies. 


Figure  5.  Initial  temperature  field  of  downstream  shoulder  as  of 
1 November  1967,  with  no  allowance  for  water-table  elevation  of  the 
tail  bay. 


Extant  reconmendatlons  elaborated  for  other  processes,  buc  ones  ap- 
proximating the  phenomenon  in  question,  may  be  utilised  for  preliminary 
calculations  [14],  [15],  [16],  [20]. 


The  area  of  the  alrswept  surface  participating  in  heat  exchange  is 
defined  as  [16]: 


where  V is  volts 

F is  area 


D ■ Z x.D  is  the  reduced  diameter 
i-1  1 1 


X is  the  amount  of  soil  having  particles  of  the 
unit  volume. 


The  volumetric  coefficient  of  heat  exchange  may  be  defined  as 


where  a is  the  local  coefficient  of  heat  exchange  between  the  air  and  the 
jointings  of  the  fill; 


is  a coefficient  allowing  for  the  difference  between  the  mean  tem- 
perature of  the  rock  fill  at  a given  point  and  the  temperature  of 
the  surface  of  a jointing. 


To  determine  the  mean  value  of  the  local  coefficient  of  heat  exchange 
a number  of  empirical  relations  are  available  [14],  [20]  by  means  of  which 
the  Nuaselt  number  may  be  determined  as  a function  of  the  Reynolds  or  the 
Grashof  and  Prandtl  number. 


With  the  Nusselt  number  known,  a may  be  defined  as 


where  X 


Coefficient  4>  , which  allows  for  the  difference  between  the  mean  tern 
perature  of  the  rock  fill  at  a given  point  and  the  surface  temperature  of  a 
jointing,  may  be  defined  as:  1 . 

d . R • 

where  Bi  is  the  Biot  number,  defined  as  Bi  - — r * 


is  the  characteristic  dimension  of  a body,  equalling  the  ratio  of 
the  volume  to  its  area. 


It  is  here  assumed  that  a regular  regime  of  thermal  interaction  between 
the  moving  air  and  the  rock  fill  arises.  The  rock  fill  represents  a two-phase 
statistical  mixture  of  jointings  of  rock  and  air.  The  following  formula  has 
been  proposed  for  determination  of  the  coefficient  of  thermal  conductivity  of 
such  mixtures  [17]: 


The  volumetric  heat  capacity  of  the  rock  fill  is  found  from  the  formula 

<^=(1  -m)0,.7„+mCv  (31) 


where  is  the  volumetric  weight  of  the  rock  fill  jointings; 

Cr  is  the  specific  heat  capacity  of  a fill  jointing. 

The  coefficient  of  temperature  conductivity  may  be  defined  as  the  de- 
rivative value 
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Relations  (28),  (29),  (30),  (31),  and  (32)  given  in  the  foregoing 
permit  determination  of  all  the  parameters  required  for  the  calculations. 
They  characterize  a rock  fill  for  which  the  thermal  regime  problem  is  solved 


7.  Method  of  Sol* 


To  determine  the  thermal  regime  of  the  downstream  shoulder  of  a rock- 
fill  dam  (simplified  model)  it  is  necessary  to  solve  equations  (12),  (22), 
and  (23)  concurrently.  Physical  quantities  6,  v,  p are  here  taken  to  be 
constant  for  the  entire  downstream  shoulder  and  correspond  to  the  mean  tem- 
perature value  of  air  in  movement. 


System  of  equations  (12),  (29),  and  (23)  can  be  solved  by  the  numerical 
method.  For  this  purpose  we  introduce  into  the  region  under  study  a rectangu- 
lar network  consisting  of  points  x ■ xn  + nAx,  y - y~  + jAy,  t ■ k*  At . 


Then  we  obtain  the  following  finite-difference  systems  for  time 


1.  For  the  equations  of  the  mean  rock  fill  temperature  (explicit 


The  temperature  of  the  points  situated  on  the  inclined  surface  of  the 
is  determined  from  the  formula 


is  the  temperature  of  the  points  adjoining  the  surface  along 


2.  For  the  temperature  of  the  air  moving  in  the  pores  (implicit  sys' 


frit if  l 


where 


k:i  - «;v  - to3j  <«■*'  - *v-‘>! 


c:r;  - <*:!  - <*w - e!- 


1 W " 1 + m\x  iy 


(KV-^->)!  + 


1,-1,  if  (KV“^+;-i>/Ay<0;  I 
1)  = 0.  if  (tv>1  - > °: 

x-1,  if  (•«.V-«-u>/Ax<0; 

y =*  0,  if 

3.  For  determination  of  the  flow  function 
(a)  in  the  Internal  nodes: 

«?  - M <«-*!. , + ,) + jv<«r;*, + *•«,) + /»,, 

i~*!-i.i) 


where 


p* 


( 


1 


4*  AjC  * Ax*  + Ay*) 


A/ 


2Ax* 


(ajc*  + Ay-) 


, A'- 


2Ay: 


1 \ 


(a^+ap) 


(36) 


(37) 


(38) 


(39) 


(b)  for  determination  of  the  flow  function  on  the  boundary  nodes 
of  the  day  (Figure  6): 


Ax 


cos  a0  + 


rj.%  + 1 Ajk+\ 

70  70 


sin  i,  - 0. 


(40) 


where  aQ  is  ths  angle  (of  inclination)  between  positive  axis  ox  and  that 
normal  to  the  day  of  the  downstream  shoulder. 
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Figure  6.  Illustrative  of  determination  of 
the  normal  derivative  on  the  day  of  the 
downstream  shoulder. 


AnedHQJf  nodcpxHQcmt 


When  the  number  of  nodes  In  the  region  In  question  is  sufficiently 
large,  it  is  virtually  Impossible  to  solve  the  algebraic  system  obtained  In 
mutual  connection.  If  In  determination  of  one  of  the  unknown  functions  for 
a given  time  use  Is  made  of  the  values  of  the  remaining  functions  from  the 
preceding  time  layer,  the  system  of  algebraic  equations  Is  found  to  be  linear 
The  concurrent  physical  processes  of  transfer  of  air  and  transfer  of  heat  by 
the  moving  air,  as  well  as  the  equation  of  propagation  of  heat,  are  In  this 
case  considered  to  be  discrete  for  each  time  step.  The  time  steps  and  the 
spatial  steps  are  selected  from  the  conditions  of  stability  and  convergence 
of  the  finite-difference  equations  obtained,  (33),  (34),  (35),  and  (38). 


Equations  (33)  and  (34)  converge  and  are  stable  if  the  following  con- 
dition Is  satisfied: 


at  all  the  points  of  the  region  under  consideration.  If  equation  (13)  Is 
approximated  by  an  explicit  finite-difference  system.  It  should  satisfy  the 
following  condition  of  stability  and  convergence  [21]  at  any  point  in  time: 


max  'V-i.y-i 


max  (%.  j—fy-i 


Equations  (38)  and  (40)  converge  and  are  stable  for  any  time  step 


The  condition  of  stability  and  convergence  of  the  system  is  thus  de- 


fined as 


maxOVy — + 


The  minimum  time  step  value  determined  from  condition  (43)  Is  sufficient 
for  the  stability  and  convergence  of  the  entire  system. 


available  to  DDU 
1 iully  legible  xeps 


If  aquation  (13)  is  approximated  by  an  Implicit  finite-difference 
system  (stable  at  any  value  At),  the  time  step  for  solving  equations  (33), 
(34),  (35),  and  (36)  must  be  selected  from  condition  (41)  only. 

Equations  (35),  (38),  and  (40)  were  solved  on  each  time  layer  by  the 
method  of  simple  Iteration. 


Using  the  distribution  of  the  temperature  of  the  air  in  the  pores  at 
• At,  we  determine  the  flow  function  for  the  following  time  step  from 


(a)  on  the  Internal  nodes 


m Vi  + n, 


Ti»,  y 


(b)  on  the  boundary  nodes  of  the  day  of  the  downstream  shoulder  (Figure 


,',♦+1+1+1 

iV/ 


The  method  of  solving  equation  (41)  is  similar  to  the  foregoing  one. 

Coefficients  A**]-,  and  were  determined  in  accordance  with  the 

n»J  «»j  ®»J  v+i  n»J 

values  obtained  for  function  7.  The  mean  temperature  of  the  rock  fill 
adopted  in  this  case  was  that  corresponding  to  time  kAT. 

The  solution  of  equation  (35)  was  assumed  to  have  been  obtained  on  a 
given  time  layer  if 

||eJi+,+,+._0»+y+l||<li  (47) 


i.a.,  if  a condition  analogous  to  (46)  was  satisfied 


The  value  of  the  air  temperature  in  the  pores  can  be  used  to  obtain 
the  temperature  of  the  rock  fill  from  (33)  and  (34). 

The  temperature  at  the  junction  of  any  two  regions  was  determined  by 
applying  the  heat  balance  equation,  allowance  being  made  for  the  difference 
in  the  thernophysical  characteristics  of  the  materials.  To  calculate  the 


thermal  regime  of  the  upstream  weighting  of  the  dan  and  base  use  was  made 
of  the  Fourier  equation  for  solids  with  constant  coefficients,  tfe  will  not 
cite  the  finite-difference  analogies  of  these  equations;  they  are  to  be 
found  in  the  corresponding  literature  [18]. 


The  calculation  work  was  accomplished  at  the  computer  center  of  the 
Krasnoyarsk  aluminum  plant,  on  an  M-20  computer.  The  assigned  program  of 
the  computer  called  for  input  of  the  initial  conditions  for  the  temperature 
of  the  rock  fill  from  punch  cards  and  referring  of  the  initial  temperature 
field  of  the  rock  fill  to  the  temperature  field  of  the  moving  air,  since 
they  were  assumed  to  be  equal  at  the  initial  time.  The  boundary  conditions 
for  the  moving  air  were  assigned  in  the  form  of  the  external  air  temperatures 
in  the  nodes  of  the  spatial  network  at  a distance  of  one  step  from  the  day. 


The  physical  parameters  were  introduced  as  a component  of  the  progr 
they  had  the  following  numerical  values: 


0.5  m; 

0.3; 

0.03  m; 

21.7  keal/h  ’ deg; 

1.293  kg/m3; 

1/273  1/deg; 

2.8  keal/m  • h * deg; 

0.005  m/h; 

■ 2.1  • 10"  3 keal/h  • deg; 

■ 1.9  keal/m  * h * deg; 

> 4.8  m; 

■ 3.0  m;  2 

« 9.81  m/sec4; 

> 0.2  (moisture  content  of  apron  and  upper  weighting) 

. 1.574  ♦ 10"6  ^ -3— C ; 


0.31  keal/m 


420  keal/m  * deg; 


External  air  temperature:  November  -20.1  °C,  December  -25°C 

January  -36.1  °C,  February  -33.4  °C.  Water  temperature  of  head  bay 


The  constant  coefficients  of  basic  finite-difference  equations  (33), 
(34),  (35),  (38),  and  (40)  were  computed  at  the  beginning  of  calculation  and 
were  stored  in  the  computer  memory.  They  were  retrieved  during  operation  by 
application  to  the  corresponding  memory  cell. 

It  took  the  longest  in  the  calculations  to  determine  the  flow  function, 
l.e.,  the  velocity  of  the  air  stream  in  the  rock  fill,  especially  at  the 
initial  time  of  calculation  and  on  change  in  the  boundary  conditions  for  air 
in  movement . 

The  temperature  of  the  rock  fill  and  the  moving  air  are  determined 
relatively  rapidly  at  each  time  step  At  - 20  hours. 

Owing  to  the  lengthy  period  of  operation,  the  possibility  of  malfunc- 
tion of  the  computer  is  not  excluded.  Hence  the  intermediate  values  of  the 
calculation  and  the  program  were  recorded  at  intervals  of  360  hours  on  a 
magnetic  drum. 

In  the  event  of  malfunction,  solution  began  with  use  of  the  most  recent 
recording.  The  results  of  the  solution  were  printed  out,  the  rock  fill  tem- 
perature at  intervals  of  180  hours  and  the  flow  function  and  moving  air  tem- 
perature at  intervals  of  360  hours. 

Figure  7 shows  the  theoretical  temperature  field  of  the  downstream 
shoulder  of  the  rockflll  dam  of  the  Vilyuy  Hydroelectric  Power  Plant.  As  is 
to  be  seen  from  the  drawing,  the  downstream  shoulder  of  the  rockflll  dam  is 
subject  to  deep  chilling. 


Figure  7.  Theoretical  temperature  field  of  downstream  shoulder  as  of 
22  February  1967  (with  no  allowance  for  water-table  elevation) . 

1.  Geometric  axis  of  dam.  2.  Apron. 


■ 


A low  temperature  with  large  gradients  Is  observed  around  the  base 
of  the  dam;  this  favors  rapid  chilling  of  the  base  of  the  dam.  Use  may 
be  made  of  this  circumstance  In  erecting  nonfiltering  rockflll  dams  on 
soft  soils.  The  elimination  of  latent  heat  during  phase  conversion  of 
the  moisture  In  the  apron  and  in  the  top  weighting,  as  well  as  the  presence 
of  a horizontal  airtight  interlayer,  exert  a warming  effect  on  the  thermal 
regime  of  the  top  of  the  downstream  shoulder  of  the  Vilyuy  rockflll  dam. 


A substantial  Influence  is  exerted  by  the  airtight  Interlayer,  which 
reduces  the  area  of  formation  of  the  air  stream  resulting  from  natural  con- 
vection. 

Figure  8 shows  the  flow  lines  of  the  air  stream  in  the  pores  of  the 
fill  on  22  February  1968.  The  concentration  of  the  flow  line  at  the  point 
of  escape  of  warm  air  demonstrates  that  the  warm  air  escapes  within  a 
limited  section  on  the  slope.  This  is  confirmed  by  the  results  of  observa- 
tions conducted  at  the  Vilyuy  dam  (Figure  3) . 
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Figure  8.  Flow  line  of  air  stream  in  downstream  shoulder  of  dam. 

1.  Geometric  axis  of  dam.  2.  Apron. 

Figures  9a,  b illustrate  the  variation  in  temperature  with  time  and 
compares  the  figures  with  the  data  of  field  observations  at  individual  points. 
The  greatest  difference  between  the  theoretical  and  the  observed  values  is 
observed  in  the  surface  zone  of  the  shoulder. 


The  reason  for  the  difference  between  the  theoretical  and  the  observed 
data  in  the  surface  zone  is  to  be  sought  in  the  influence  of  the  variations 
in  external  air  temperature,  the  Initial  conditions  of  the  problem,  the 
heterogeneity  of  the  rock  fill,  and  the  Influence  of  the  heat  insulating 
and  airtight  properties  of  snow  on  the  downstream  slope  of  the  dam. 
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Figure  9.  Variation  In  temperature  at  various  points  of  the  downstream 
shoulder  with  time,  (a)  at  point  N°  5;  (b)  at  points  N°  14  and  18; 

1.  theoretical  curve  without  allowance  for  water-table  elevation  of  tail 
bay.  2.  same  with  allowance  made  for  water-table  elevation  of  tail  bay. 

3.  according  to  data  of  field  observations.  4.  November  1967.  5.  Decem- 
ber. 6.  January.  7.  February  1968.  8.  hours.  9.  thermistor. 

9.  Thermal  Regimes  of  the  Downstream  Shoulder  of  a Rockfill 
Dam  with  a Day  Having  an  Airtight  Surface  Cover 

Such  problems  can  arise  if  the  surface  of  the  downstream  shoulder 
is  covered  with  a thick  layer  of  snow  or  other  airtight  materials.  The 
method  of  solving  the  problem  remains  unchanged  in  this  case;  only  the 
boundary  conditions  of  the  problem  undergo  variation. 

In  solving  the  transfer  function  it  is  necessary  to  assume  the  flow 
function  to  equal  zero  over  the  entire  outline  of  the  downstream  shoulder. 

The  temperature  of  the  day  of  the  downstream  shoulder  is  determined 
from  the  equation  ^ 

} ~dn^  “ef^ext-  (54) 

where  X is  the  coefficient  of  thermal  conductivity  of  the  rock  fill 
without  convection; 


1 

1 + 1 

a X. 


is  the  effective  coefficient  of  heat  exchange  of 


the  rock  fill  with  air  [13]; 

a is  the  coefficient  of  heat  exchange  of  the  fill  with  air; 
l is  the  airtightness  thickness  of  the  material; 

X is  the  coefficient  of  thermal  conductivity  of  the  airtight  material. 


22 


The  temperature  of  the  pore  air  Is  everywhere  determined  by  finite- 
difference  equation  (35).  The  temperature  of  the  air  in  the  pores  of  the 
day  of  the  downstream  shoulder  is  to  be  adopted  in  this  instance  as  the 
origin  of  the  temperature  level. 

The  numerical  values  of  the  parameters  employed  in  the  calculations 
were  assumed  to  be  the  same  as  those  previously  used,  with  the  exception 
of  the  coefficient  of  permeability  (k'  *0.34  • 10“  ^ m2)  and  tjje  time  step 
(At  ■ 10  hours). 

The  thermal  regime  of  the  base  of  the  dam  was  not  considered.  The 
boundary  temperature  in  the  base  was  assumed  to  equal  zero  degrees.  The 
initial  condition  of  the  problem  is  also  slightly  different  (Figure  10). 
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Figure  10.  Initial  temperature  field  of  downstream  shoulder. 
1.  geometric  axis  of  dam. 
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Figure  11.  Theoretical  temperature  field  of  downstream  shoulder  with 
airtight  cover.  1.  First  layer  of  filter.  2.  Second  layer  of  filter 

Figure  11  illustrates  the  theoretical  temperature  field  of  the 
downstream  shoulder  when  an  airtight  cover  is  present,  and  Figure  12  the 
temperature  field  problem  solved  with  no  airtight  cover. 
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Figure  12.  Theoretical  temperature  £ield  of  downstream  shoulder  without 
airtight  cover.  1.  geometric  axis  of  dam.  2.  Apron 


Comparison  of  the  two  temperature  fields  (Figures  11  and  12)  referred 
to  the  same  point  in  time  demonstrates  that  the  temperature  of  the  downstream 
shoulder  is  higher  when  an  airtight  cover  is  present. 

The  presence  of  an  airtight  cover  reduces  the  speed  of  movement  of  the 
pore  air  by  a factor  of  2 to  3.  The  mean  temperature  of  the  downstream  shoulder 
varies  at  a slower  rate  when  an  airtight  cover  ia  present  than  when  one  is 
absent  (Figures  11  and  12). 

The  airtight  cover  of  the  entire  slope  of  the  downstream  shoulder  can 
be  utilized  for  thermal  regulation  of  rockfill  dams.  The  presence  of  an  air 
insulating  cover  during  the  winter  period  alone  fosters  general  elevation  of 
the  temperature,  and  during  the  summer  period  lowering  of  the  temperature,  of 
the  downstream  shoulder  of  a dam. 

Separated  closed  air  streams  are  formed  in  the  downstream  shoulder  of 
a rockfill  dam.  Convective  heat  exchange  between  individual  streams  is  ab- 
sent, and  exchange  of  heat  between  them  is  effected  only  by  way  of  thermal 
conductivity. 

In  conclusion  it  is  to  be  noted  that  we  have  not  made  a full  study  of 
the  influence  of  an  airtight  covering  of  the  day  of  a downstream  shoulder  on 
the  thermal  regime. 

The  problem  must  be  solved  for  a more  prolonged  period  (2-3)  years  in 
order  to  obtain  more  reliable  results. 


Conclusions 


1.  A study  has  been  made  of  the  thermal  regime  of  the  downstream 
shoulder  of  the  rockfill  dam  of  the  Vilyuy  Hydroelectric  Power  Plant  on 
the  basis  of  design  relationships  (12),  (22),  and  (23). 


2.  The  results  of  cslculation  and  the  data  of  field  observations 
have  been  found  to  be  in  close  agreement. 


3 Certain  differences  between  them  are  to  be  ascribed  to  determine’ 
tlon  by  way  of  approximation  of  the  parameters  entering  into  the  design  re- 
lations (porosity  of  the  fill,  coefficient  of  permeability  of  the  rock  fill 
diameter  of  the  jointings,  etc). 


4.  The  results  of  the  study  can  be  employed  in  the  design  of  rock- 
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